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Recently a four-loop calculation of the gauge /3-function (f3 g ) in an N — 1 super- 
symmetric gauge theory was presented [[I]. The computation was performed within the 
usual supersymmetric dimensional regularisation scheme (DRED), which is standard for 
perturbative calculations in supersymmetry and which consequently would be the most 
convenient for any future phenomenological applications. The result was derived some¬ 
what indirectly, by starting from a partial calculation for the abelian case. The coupling 
constant redefinition relating the abelian result to the corresponding exact NSVZ result [2| 
could then be constructed, and this redefinition was then extended to the non-abelian case 
by imposing the vanishing of the /3-function beyond one loop in the case of AT — 2 super- 
symmetry. However, one possible redefinition which was present in general, but vanished 
for the abelian case, could not be determined, and consequently the four-loop result /?3 
contained one undetermined constant. Our purpose here is to predict the value of this 
constant by using improved Pade approximation techniques which have recently been de¬ 
veloped and applied to the case of standard QCD0. In that case a prediction for (3^, was 
found to be in good agreement with a subsequent analytic calculation]^], especially when 
the contributions involving quartic Casimir group theory structures (which arise for the 
first time at four loops) are omitted in making the comparison. In fact these structures 
simply do not arise in the supersymmetric case which concerns us here. This observation, 
allied with the fact that we have more information about (3% than was available to Ellis et 
al in the QCD case, makes us optimistic regarding the accuracy of our prediction. 

The method is as follows: For a perturbative series P(x) = S n x n , the Pade 

approximant P[n/m]( x ) is a ratio of polynomials A^{x) and of degree N and M 

respectively, chosen so that 

p [K /Mi = Tdh = p(i) + o(x A ' +M+I ). (i) 
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It can be argued that in the case of QCD (and presumably also of supersymmetric QCD) 
the relative error 
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has the asymptotic form 
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as N —> oo, for fixed M, where A is a constant and L — N + M + aM + b. We will be 
concerned with [0,1] and [1,1] Pades, so M = 1; and we chooseQ a, b so that a + b — 0. 
The Asymptotic-Pade Approximant Prediction (APAP) is then given by: 
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If Sq, Si and 57 are known, then the prediction for S 3 is obtained as follows. Firstly by 
matching to P(x) up to linear terms as in Eq. (|l[), we obtain 


P [o,i] = 
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giving as a prediction for the coefficient of the x 2 term £' Pade = -<A. Then from Eq. (|2|) we 
hnd §2 — g o g 2 — 1. From Eq. (§) we deduce A = —62 , and thence ^3 = —\A = ^ 62 - The 
[1,1] Pade prediction for S 3 is easily derived to be S 3 ade = |£, and we then use Eq. 
give an improved estimate S^ PAP incorporating ^ 3 . 

The gauge /5-function is given by 


to 


OO 

Pg = 9^2PnX n+1 , ( 6 ) 

n =0 
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where x — and (3 n corresponds to the (n + l)-loop result in perturbation theory. 

The procedure given above yields a prediction for the four-loop contribution /S 3 . The 
final feature of the method is to perform the above process for a range of values of Nj-. the 
number of flavours, and then to match the results to a cubic polynomial in Nf. In the QCD 
case this enabled the authors of Ref. |§ to incorporate as an extra piece of information 
the known coefficient of Nj, obtained by large- Nf calculations. In our present case we 
shall be able to exploit the fact that we know the form of P 3 as a function of Nf up to 
a single parameter. The results of Ref. [lj were presented for a general theory; for the 
present purposes, however, we specialise to supersymmetric QCD, obtaining [Q: 
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where N c is the number of colours, and 
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Here k — 6£(3) and a is the constant which we hope to determineB. 

Note that T> does not involve a. It is the leading-iVj contribution at four loops and 
could, of course, have been obtained from the large-iVy results presented in Refs. |[|. |7]|. 
We may therefore proceed in much the same way as did the authors of Ref. |Q. We first 
compute S^ FAF as described above, and then after determining A. B. C from the fit to 
Eq. ([ 7^) we can obtain three distinct “predictions” for a, a (A), a(B) and a(C). Naturally 
a test of the approach is the extent to which these predictions agree and are insensitive 
to the input value of N c . As in Ref. we use input values Nf — 0,1, 2, 3, 4, and used a 
value for A, (A), obtained by averaging A(Nf ) for the input values; the outcome is not 
very sensitive, in fact, to whether we use (A) or A(Nf). The results of this procedure are 
shown in Table 1: 


N c 

a{A) 

a(B) 

«(C) 

ct(C) (k — 0) 

2 

2.56 

2.34 

-0.31 

1.62 

3 

2.43 

2.42 

0.50 

2.38 

4 

2.46 

2.47 

0.63 

2.48 

5 

2.45 

2.46 

0.74 

2.57 


Table 1: APAP predictions for the unknown parameter a. 

The results for a(A),a(B) are remarkable for their consistency and stability. In the 
fourth column of the table we show the result for a(C) consequent on omitting from /R 

1 It is also possible to apply the APAP method to the cases of QED or supersymmetric QED, for 
which the complete four-loop results can be extracted from Refs. || or JTJ respectively. However, 
the form of the series in these cases is not conducive to the accuracy of the approximation and 
the four-loop predictions are somewhat less impressive. We hope to discuss these and other 
applications of the APAP method elsewhere. 
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the terms proportional to k. It could be argued that this is natural since, given that £(3) 
occurs for the first time at four loops, these terms cannot be accurately produced by the 
the Pade. (This is similar to the argument made with respect to the quartic Casimir terms 
in Ref. ||.) In any case, our results clearly suggest a value for a of around 2.45. 

An alternative approach is as follows. The unknown parameter a in Eq. (|8|) is one of 
three parameters that were introduced in Ref. [U, each being the coefficient of a certain 
coupling constant redefinition. Two of these parameters were in fact determined in Ref. [![]; 
but we can test the reliability of the APAP method by reintroducing one of them and 
comparing the APAP prediction with a f3s now dependent on two parameters. Thus we 
replace /?3 by /R where: 


ft -> ft = A + BN, + CNj + VN'j 


(9) 


and 
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In the notation of Ref. |l[], a = ot\ and 8 = — | — ot 2 ■ In Ref. [|l| it was shown that a 2 = — 
and so we hope to find that 5 = 0. For fixed a we determine three distinct “predictions” 
for 8 by fitting to a polynomial as before. As an example, for Ay. = 4 we obtain: 


8 (B) = 7.85cr - 19.24 
5(C) = 19.38a-48.16 
8 (V) = 51.94a- 126.87. 


(11a) 

(116) 

(11c) 


If we omit k terms then Eqs. Qllrj) , ( |llc|) are unaffected, but Eq.( |l 1 /{ ) becomes 
5(C) = 19.38a-47.41. 

When we plot 5(B), 8(C) and 8 (V) against a, as in Fig. 1, the results are quite striking. 
We see that 8 (B), 8(C) and 8 (T>) are zero for very similar values of a, around 2.44-2.49. 
So the predicted values for 8 (B), 5(C) and 8 (T>) agree rather well with the correct value 
5 — 0 if 2.44 < a < 2.49, and this may be interpreted as a prediction for a. 
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The convergence is again even better if we omit the term in k. as depicted in Fig. 2, 
where each 5 crosses the axis at a — 2.44. As N c is increased, we find that the convergence 
of the 5s improves in the case where we retain k; in fact for large N c all the 5s are zero for 
a — 2.43 irrespective of whether or not we retain k. 

What is our final prediction for o? The relatively simple relationship between the 
NSVZ /3-function and the corresponding DRED one as explored in Ref. JT|] suggests that 
a is a simple fraction; given our results here, then we would expect a — 12/5 or perhaps 
a — 5/2. (a = 17/7 is even closer, but experience suggests that this number is unlikely to 
emerge from a perturbative calculation.) The result a = 2£(3) is also possible; in which 
case the apparent better convergence produced by excising k from the comparison would 
need to be dismissed as coincidental. 
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